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We compare the theoretical predictions of the previous article on the field dependence of the mag- 
netic spectrum in anisotropic two-dimensional and Dzyaloshinskii-Moriya layered antiferromagnets 
[L. Benfatto and M. B. Silva Neto, cond-mat/0602419 , with Raman spectroscopy experiments in 
Sr2Cu02Cl2 and untwinned La2Cu04 single crystals. We start by discussing the crystal structure 
and constructing the magnetic point group for the magnetically ordered phase of the two compounds, 
Sr2Cu02Cl2 and La2Cu04- We find that the magnetic point group in the ordered phase is the mmm 
orthorhombic group, in both cases. Furthermore, we classify all the Raman active one-magnon ex- 
citations according to the irreducible co-representations for the associated magnetic point group. 
We find that the in-plane (or Dzyaloshinskii-Moriya) mode belongs to the DA g co-representation 
while the out-of-plane (XY) mode belongs to the DB g co-representation. We then measure and 
fully characterize the evolution of the one-magnon Raman energies and intensities for low and 
moderate magnetic fields along the three crystallographic directions. In the case of La2Cu04, a 
weak-ferromagnetic transition is observed for a magnetic field perpendicular to the Cu02 layers. 
We demonstrate that from the jump of the Dzyaloshinskii-Moriya gap at the critical magnetic field 
H c ~ 6.6 T one can determine the value of the interlayer coupling J±/J — 3.2 x 10 _J , in good 
agreement with previous estimates. We furthermore determine the components of the anisotropic 
gyromagnetic tensor as = 2.0, g b s = 2.08, and the upper bound = 2.65, also in very good agree- 
ment with earlier estimates from magnetic susceptibility. For the case of S^CuChCh, we compare 
the Raman data obtained in an in-plane magnetic field with previous magnon-gap measurements 
done by ESR. Using the very low magnon gap estimated by ESR (~ 0.05 meV), the data for the 
one-magnon Raman energies agree reasonably well with the theoretical predictions for the case of a 
transverse field (only hardening of the gap). On the other hand, an independent fit of the Raman 
data provides an estimate for g s ~ 1.98 and gives a value for the in-plane gap larger than the 
one measured by ESR. Finally, because of the absence of the Dzyaloshinskii-Moriya interaction in 
Sr2Cu02Cl2, no field-induced modes are observed for magnetic fields parallel to the Cu02 layers in 
the Raman geometries used, in contrast to the situation in La2Cu04. 

PACS numbers: 74.25.Ha, 75.10.Jm, 75.30. Cr 



In the preceding article^ we investigated the field dependence of the magnetic spectrum in anisotropic two- 
dimensional and Dzyaloshinskii-Moriya layered antiferromagnets. The first case is relevant to the understanding of 
the magnetic properties of Sr 2 Cu02Cl2, which is a body-centered tetragonal S = 1/2 antiferromagnets with IA/mmm 
structure 3 - and D4/1 point group in the paramagnetic phase, T > Tm- Because of its tetragonal character, at the clas- 
sical level there should be no in-plane anisotropies present in S^CuCbC^. Furthermore, because of its body-centered 
structure, J± is strongly frustrated^ However, such perfect frustration can be removed by quantum fluctuations due 
to the spin-orbit interaction, and indeed a small in-plane anisotropy is present in S^CuC^C^^ determining a spin 
easy-axis at low temperatures, and giving rise to a very small in-plane spin gap- 6 ' For this reason, the magnetism in 
Sr2Cu02Cl2 can be fairly well described by the following two-dimensional square-lattice spin-Hamiltonian 
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In the above expression, S^' b ' c are the components along the crystallographic axes a, b, c (see Fig. of the Cu ++ 
spin, J is the planar antiferromagnetic superexchange, and a a and a c are, respectively, parameters that control the 
in-plane and XY anisotropics. It should be emphasized here, once more, that at the classical level a a — because 
the crystal structure is body centered tetragonal, 7 A nonzero a a can however be effectively obtained once quantum 
fluctuations (mostly from the spin-orbit coupling), which lift the frustration, are taken into account by considering 
the mean field effect of the neighboring layers^ 




PIG. 1: Tetragonal magnetic unit cell of SnCuChCh for T < Tn (only Cu ++ , 0~~ and CP ions are shown for clarity). For 
the sake of unifying notations between S^CuChCh and La2Cu04, in what follows we shall use the (abc) coordinate system 
above where a \\ x' , b || y' , and c || z. Since the crystal structure is tetragonal, we have a = b < c. The Cu ++ spins are confined 
to the (ab) plane (CuO-2 layers) and are oriented along the (110) direction in the (xyz) coordinate system, or, equivalently, 
parallel to the b axis. In the paramagnetic phase, the point group is the IA/mmm tetragonal group, but, as we shall discuss 
soon, below Tn the symmetry is lowered because of the antiferromagnetic ordering of the spins. 



Conversely, any realistic model for the magnetism of La2Cu04 that takes into account the tilting of the oxygen 
octahedra should incorporate both a Dzyaloshinskii-Moriya (DM) interaction and also the interlayer coupling J± ^ 0£ 
In the low-temperature orthorhombic phase of La2Cu04, T < 530 K, the crystal has the Bmab structure with the D2h 
point group in the paramagnetic phase, T > Tn- The full Hamiltonian that incorporates the Dzyaloshinskii-Moriya 
and XY interactions allowed by symmetry, as well as the interlayer coupling, reads 

h = j± s '" • s ™ +1 + Hsi [ s?n ' Dm ] ' ( 2 ) 

m m 

where S m and D m represent, respectively, the Cu ++ spins and Dzyaloshinskii-Moriya vectors at a generic lattice 
position of the mth layer. The sum is over the Hamiltonian for a single layer 

H sl [S, D] = J J2 S * ' S j + Yl D ^ ' & x S J') + J2 S * ' ^ ^ ' S J ' ( 3 ) 

where Dy and T ij are, respectively, the DM and XY anisotropic interaction terms that arise due to the spin- 
orbit coupling and direct-exchange in the low-temperature orthorhombic (LTO) phase of La2Cu04 (see Fig. Eland the 
preceding article for a proper definition of these quantities) . It should be noted here that due to the peculiar staggered 
pattern of the tilting angle of the oxygen octahedra in neighboring layers, the Dzyaloshinskii-Moriya vector alternates 
in sign from one layer to the other, D™ B Ac = (— \ ) m T) ab ,ac ■ Moreover, since the unit cell is body centered, the 
coupling Jj. in Eq. 1J2J connects the spin at position (0,0,0) to the one at (1/2,0,1/2) in the LTO coordinate system 
(see also Fig.EJ. 

As it has been discussed extensively in the preceeding article^, the DM interaction leads to a quite unconventional 
response of the system to an external magnetic field. Indeed, due to the DM interaction the spins develop small out- 
of-plane ferromagnetic moments along with the staggered ones characteristic of the AF order, which in turn couple 
to the external field leading for example to unusual magnetic susceptibility anisotropies, as measured at small fields 
in La2Cu04iS*iSiiiid At larger field values one can observe a ferromagnetic ordering of these moments along the c 
axis for magnetic fields applied perpendicular to the CuCh layers ^i 1 ^ 11 or to a two-step spin-flop of the staggered 
moments for an applied in-plane fielcUfiiii. This physical picture, which was discussed in the previous work of Refs. 

has been confirmed by the calculations reported in Ref. pj, where additional results concerning the (H,T) 
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phase diagram and Raman response at finite magnetic field have been discussed. This same approach turned out to 
be quite successfull recently in demostrating that the Dzyaloshinskii-Moriya interaction is behind the appearance of a 
field-induced mode in the one-magnon Raman spectrum in La2CuC>4 for longitudinal magnetic fields, as a consequence 
of a rotation of the spin-quantization basis which was first suggested in Ref. fill theoretically explained in Ref. 0, 
and quite recently observed in neutron diffraction in Ref. 0. Thus, the purpose of the present paper is to directly 
compare the one-magnon spectrum measured in both La2Cu04 and S^CuChC^ with the calculations presented in 
Ref. Q, as far as both the position and intensity of the Raman peak is concerned. As we shall see, the excellent 
agreement between the thoery and the experiments allows from one side to establish on firmer grounds the general 
picture of the anomalous effect of the DM interaction in the La2Cu04 system, and from the other side to estimate 
some relevant physical quantities as the spin gyromagnetic ratio. 




PIG. 2: (Color online) - Orthorhombic magnetic unit cell of La2Cu04 for T < Tjv (only Cu ++ and 0~~ ions -except the 
apical ones- are shown for clarity). We use the (abc) orthorhombic coordinate system, where a ^ b c. The Cu ++ ions order 
antiferromagnetically parallel to the b orthorhombic axis but are also canted out of each CuCh layer, and are thus confined to 
the (6c) plane. Observe that the canting pattern is staggered along the c orthorhombic direction. 




FIG. 3: Left: filled (small) circles represent the Cu ++ ions. The hatched larger circles represent the oxygen O ions that are 
tilted above the Cu02 layer while the open larger circles are O ions that are tilted below the Cu02 layer, see also Fig. E] 
Right: the staggered pattern of the Dzyaloshinskii-Moriya vectors, represented by open arrows along the Cu-Cu bonds, follows 
from the staggered pattern of the tilting of the oxygen octahedra. Small arrow: in-plane projection of the Cu ++ spins in the 
ordered phase. 
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II, MAGNETIC GROUP ANALYSIS FOR La 2 Cu0 4 AND Sr 2 Cu0 2 Cl 2 

Before introducing the experimental setup, a preliminary discussion is needed about the point group of La2Cu04 
and Sr 2 Cu02Cl2 both in the paramagnetic and antiferromagnetic phases, in order to establish the irreducible co- 
representations of the Raman tensors. While the paramagnetic case has been already extensively discussed in the 
literature (see for example Ref. ^3 an d references therein) , a discussion of the magnetic phase is still missing, and we 
shall briefly review here the basic steps needed to make this analysis. 

Let us consider a crystal with a certain symmetry group, G, and its set of allowed symmetry group operations in 
the paramagnetic phase, T > Tjy. When the system orders antiferromagnetically, T < Tjy, usually the symmetry is 
lowered because now not all the sites are equivalent, and the corresponding magnetic group must be determined^ In 
general, two situations can occur: 

(a) Let us call H the restricted unitary subgroup of G containing the symmetry operations still allowed below Tjy. 
If H is a subgroup of G of index 2 (i.e. H contains exactly n/2 elements of G) then all the remaining elements 
of G — H can nevertheless be promoted to allowed symmetry operations when combined with the time-reversal 
operation 9, and are thus called anti-unitary elements. One can then identify, for T < Tjsr, the magnetic point 
group corresponding to the classical point group G as 

M = H + 0(G-H). (4) 

(b) However, it is possible that below Tjv the unitary operations still allowed form a subgroup 7i of index larger than 
2 for G, but corresponding to a subgroup of index 2 of a different classical group Q (i.e. H contains r/2 < n/2 
elements, where r is the dimension of the group Q). In this case the magnetic group is identified by H and the 
anti-unitary group 9(Q — 7i), so that 

M = H + e(G -H). (5) 

A. Magnetic point group for La 2 Cu04 

The crystal structure of the low temperature orthorhombic phase of La2Cu04 is the Bmab structure in Fig.[2j which 
has as unitary group the Din point group in the paramagnetic phase T > Tjv. Above the Neel ordering temperature 
the allowed symmetry operations of the Dih group are 

-D2/1 = 1, 1, 2 a , 2b, 2 C , 2 a , 2b, 2 C , (6) 

where the 8 elements have their usual meaning^: 1= identity; 1= inversion through the symmetry center, which is 
in this case the central Cu ++ ion in position 2 of Fig. [21 so that (a, b, c) — > (—a, —b, — c); n z = rotation of 2ir/n around 
the z axis; n z = inversion through the symmetry center followed by a 27r/n rotation around the z axis (note that 2 Z 
corresponds also to a reflexion with a mirror perpendicular to the z axis). 

Although it may seem at first that neither 2b nor 2 C are symmetry operations, in both cases the final atomic 
configurations can be brought back to the original one by a translation of half of the diagonal along the (Oil) 
direction (in the (abc) coordinate system). In this sense, the octahedra in position 1 (front-bottom- left corner) is 
brought to position 2 (central ion) and the one from position 2 is brought to position 3 (back-top- right corner). Such 
translation is allowed because there is no way to distinguish the corner and central Cu ++ ions in the crystal. 

Below Tjy, the Cu ++ ions order antiferromagnetically in the pattern shown in Fig. We can verify that the 
remaining allowed unitary symmetry operations are 

H = 1, 1, 2b, 2b, 

bearing in mind that a half translation along the (Oil) diagonal is allowed. We immediately conclude that La2Cu04 
belongs to the case (a) discussed above, where the subgroup of allowed unitary operations is of index 2, see Eq. J3J, 
We can now construct the magnetic group of the ordered phase of La 2 Cu04 by combining all the G H unitary 
operations that are not allowed below T/v with the time reversal operation 9 (represented in what follows by an 
underline), thus forming anti-unitary operations. The allowed unitary + anti-unitary operations for La2CuC>4 below 
Tjv are 



M = H + 6(G - H) = 1, 1, 2 a , 2 b , 2 C , 2 a , 2 b , 2 C , 



(7) 
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such that the magnetic point group for La2Cu04 below Tjv is the 

M = mmm 

orthorhombic group, which also has 8 elements. The Raman tensors for such magnetic group are given in terms of 
the co-representations (see Cracknel!^) 

/ A LB \ / C* \ 

DA g = I iD E , DB g = if , (8) 
\ / / \G iH J 

where A, B, C, D, E, F, G, H, I are unconstrained real numbers. 



B. Magnetic point group for Sr2Cu02Cl2 

The crystal structure of the high temperature tetragonal (HTT) phase of S^CuC^C^ is the I4/mmm structure in 
Fig.QJ which has as unitary group the D^h point group in the paramagnetic phase T > Tjy. Above the Neel ordering 
temperature the 16 allowed symmetry operations of the D^h group arei& 

1> I5 2 X 7 ^y, 2^, 2^' j 2yl j 2^, 2y, 2^, ^x'l 2yf , it 4_g , it 4^; , 

where we used the reference system of Fig. ^ As we can see, the tetragonal character of the unit cell allows for a 
4-fold axis, the c-axis, which is perpendicular to the CuC>2 planes (see Fig.QJ. 

In the antiferromagnetic phase the spin easy-axis is along the y' or (010) direction in the (abc) coordinate system, 
see Fig. n] Thus, one can easily verify that, below T/v, 2 x , 2 y , and ±4 2 are no longer symmetry operations, not even 
when supplemented by the time reversal operation 9. The only allowed unitary operations in the Neel ordered phase 
of Sr 2 Cu02Cl2 are 

/T 1, 1, lyl , lyl , 

again bearing in mind that a half translation along the (Oil) diagonal is allowed. We immediately conclude that 
Sr 2 Cu02Cl2 belongs to the case (b) discussed above, as in Eq. (0, where the subgroup H of allowed unitary 
operations is a subgroup of D^h with index larger than 2, but it is a subgroup of Q — Z?2h (mmm) with index 2. 
We can now construct the magnetic group of the ordered phase of S^CuC^C^ by combining all the Q — H unitary 
operations that are not allowed below Tjv with the time reversal operation 9, thus forming anti-unitary operations. 
The allowed unitary + anti-unitary operations for Sr2Cu02Cl2 below Tjv are 

M = H + 6(Q -H) = 1, 1, 2 X ,, 2 y ,, 2 Z , 2 X „ 2 y ,, 2 Z , (9) 

and thus the magnetic group, with only 8 elements, is again the orthorhombic mmm group, just like the case of 
La2Cu04i2£ In fact, when expressed in terms of the (abc) coordinate system of Fig. such that a || x 1 , b \\ y' , and 
c || z, Eq. ® can be written: 

M = 1, 1, 2 a , 2 b , 2 C , 2 a , 2 b , 2 C , 

which coincides with Eq. above. Thus, to unify the notation we shall use in the following the (abc) coordinate 
system for both La2CuC>4 and Sr2Cu02Cl2 while discussing the antiferromagnetic phase. Moreover, since the magnetic 
group is the same, the Raman tensors are given by Eq. ||HJ for both systems. 



III. INELASTIC LIGHT SCATTERING BY MAGNONS IN La 2 Cu0 4 AND Sr 2 Cu0 2 Cl 2 

One of the possible mechanisms for the inelastic scattering of light by magnetic excitations in crystals is an indirect 
electric-dipole (ED) coupling via the spin-orbit interaction*^ Such a mechanism has been in fact used to determine the 
spectrum of magnetic excitations in many different condensed matter systems like fluorides, XF2, where X is Mn 2+ , 
Fe 2+ , or Co 2+ , 21 inorganic spin-Peirls compounds, CuGeOar^ and the parent compounds of the high-temperature 
superconductors^ The Hamiltonian representing the interaction of light with magnons can be written quite generally 



^=^E^(r)E J; 

r 
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where Eg and Ej are the electric fields of the scattered and incident radiation, respectively (a T is the transposed of 
the a vector) and x( r ) is the spin dependent susceptibility tensor. We can expand x( r ) in powers of the spin-operators, 
S, as 

X aP (r) = X f(r) + J2 K a ^{v) + ]T G a ^{r)S v {v) + . . . , 

where fi, v = x, y, z label the spin components. The lowest order term Xo /3 (r) is just the susceptibility in the absence 
of any magnetic excitation (it corresponds to elastic scattering), and it will be neglected in what follows. The second 
and third terms can give rise to one-magnon excitations because they can be written as S' ± (r) and 5' z (r)S' ± (r), 
respectively, where here z is the direction of the spin easy-axis. The intensity of the scattering, as well as the selection 
rules, will be determined by the structure and symmetry properties of the complex tensors K and G. 

The ED Hamiltonian that describes the k = one-magnon absorption/emission on sublattice A/B can now be 
written in terms of the x,y,z components of the sublattice magnetization, Mi = (Sj A — Si B )/2, as 

a ed = £ {e£ x x e 7 M f + e£ x ^e 7 m\ } , (io) 

i 

where the matrices \ x an d X v are both written in terms of the original K and G tensors. Here we made the usual 
mean-field assumption (S* A ) = —(S? g ) = — S and we dropped terms of the type S^' y + S^ v since these do not 
contribute for the k = scattering. 



A. Magnetic excitations and selection rules for La2Cu04 and Sr2Cu02Cl2 

It is very important to emphasize here that, differently from two-magnon Raman scattering, where the Raman 
response does not depend on the direction of the spin easy-axis, the one-magnon Raman response does. Indeed, it 
has been shown by Fleury and Loudon2i that one of the incoming/outgoing components of the electric field must 
always lie in the direction of the easy axis, and the other one in the perpendicular plane, oriented in the direction 
of the mode that one wants to probe. In the specific case of La2Cu04 and S^CuC^C^ that we are considering, the 
easy axis is along b. Thus, the only nonvanishing matrix elements for \ x = x a , which corresponds to the in-plane 
or a (DM) mode, are those that mix the spin easy-axis, b, with a, while for x y = X c > which corresponds to the 
out-of-plane or c (XY) mode, are those that mix b and c. This information allows us to establish: (i) to which 
magnetic-group co-representation the one-magnon modes belong, and (ii) the precise scattering geometries needed to 
observe the associated one-magnon Raman modes. For example, it tells us that, for either S^CuC^C^ or La2CuC>4, 
the out-of-plane mode (also referred to as XY mod o 6 i 15 ) belongs to the DB g co-representation of the magnetic point 
group mmm 



r = X XY = | LP J , (11) 
iff 

with C = G = 0. On the other hand, the in-plane mode (also referred to as DM mode for L^CuO^A***^), belongs to 
the DA g co-representation of the magnetic point group mmm 



DM 



iB 

LD | , (12) 




with A = E = I = 0. The precise numerical evaluation of the remaining elements B,D,F,H requires a more detailed 
microscopic calculation of the electric-dipole induced transitions and spin-orbit coupling in second order perturbation 
theory, for each system, and this is beyond the scope of this article. 



B. One-magnon Raman response 

According to the electric-dipole Hamiltonian ltTfl)l . one-magnon Raman scattering probes the long- wavelength spin 
excitations Mi around the staggered spin configuration realized in the antiferromagnetic state. These correspond to 
the (k = 0) value of the usual spin-wave dispersions, where k is measured with respect to the Q = (w, ir) vector of the 
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antiferromagentic ordering. In an isotropic Heisenberg antiferromagnet the spin modes are soft at long wavelength, 

cjk — ck, but for anisotropic spin-spin interactions a gap appears at k = 0, cjk = J uj^ c + c 2 k 2 , where uj a . c are related 

to the anisotropy parameters (a a ,a c in Eq. Q, or to D + ,T in Eq. QJ). Starting from the Hamiltonian IjlOJl . the 
Raman response can be obtained using Fermi's golden rule and, for Stokes scattering, we havei& 

x(u) = [n fl (w) + 1] ||n Q | 2 A(o,^) + |n c | 2 A(o,^)} , (13) 

where T(w) is the Raman intensity, n B (uj) = (e^ — is the Bose function and *4 a ,c(0, w) is the spectral function 
of the k = transverse a/c components of the antiferromagnetic (staggered) order parameter. The projectors 

n a , c = e£ x 0,c Ei (14) 

are given in terms of the Raman tensors % a ' c presented above. The properties of A a , c have been extensively discussed 
in the preceeding paper i In particular, it has been shown that the spectral function A a , c (0,oj) is peaked at the energy 
u a , c {H) of the magnon gaps in magnetic field, with an intensity I a , c which also depends on H. For La2Cu04, uo a 
corresponds to the Dzyaloshinskii-Moriya gap, while for S^CuC^C^ it corresponds to the in-plane gap, as discussed 
in the introduction. For both compounds, lj c corresponds to the XY anisotropy gap. 



IV. EXPERIMENTAL SETUP 



Single crystals of La2Cu04 and S^CuChCb 24 were measured in a backscattering geometry with the incoming 
photons propagating along the c crystallographic axis. The polarization configurations are denoted by (E/,Es) with 
Ej/s representing the direction of the incoming/scattered electric field. The crystals were mounted in a continuous 
flow optical cryostat and the Raman spectra were taken using about 5mW power and the wavelength A = 647.1 nm 
from a Kr + laser. The measurements in magnetic fields were taken with the cryostat inserted in a room temperature 
horizontal bore of a superconducting magnet. The orthorhombic axes of the La2Cu04 sample were identified by X-ray 
diffraction. The data from S^CuCbC^ crystal were taken from a freshly cleaved (ab) surface but in this case the 
in-plane axes directions were not determined. 




FIG. 4: (Color online) One-magnon Raman response for La2Cu04 at T = 10 K in the (RL) circular polarization configuration, 
where only the Dzyaloshinskii-Moriya gap is Raman active. In (a) the gap only hardens, in (6) the gap only softens, while in 
(c) it first softens and then hardens above the weak-ferromagnetic transition at H ~ 6.6 T (see discussion in the text). 



We first present the one-magnon Raman response in La2Cu04 for the (RL) polarization configuration. In this 
geometry, the electric-field of the incident light is circularly polarized rotating clockwise, ~Ef = (x a — ijc&)/v2, while 
the electric-field of the scattered light is circularly polarized and rotating anti-clockwise, Eg = (x a + ix&)/V2. Here 
x a and X(, are unit vectors along the a and b directions respectively. According to our earlier discussion on the 
classification of the magnetic excitations in La2Cu04, the (RL) in-plane polarization configuration is the adequate 
Raman geometry in order to observe the Dzyaloshinskii-Moriya gap, because it probes directly the nonvanishing 
element B,D of the \ DM matrix 02}. The results are presented in Fig.Q] In (a) the magnetic field is applied along 
the a axis (transverse field). We observed a monotonic hardening of the gap with increasing magnetic field. In (b) 
the magnetic field is applied along the b axis (longitudinal field) . We observed a monotonic softening of the gap with 
increasing magnetic field. Finally, in (c) the magnetic field is applied perpendicular to the CuC>2 layers (transverse 
field). We observed first a softening of the gap with increasing field, a jump at a field of approximately 6.6 T, and 
finally a hardening with increasing field. 
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FIG. 5: (Color online) One-magnon Raman response for La2Cu04 at T — 10 K in the (RR) circular polarization configuration 
for a magnetic field applied along the orthorhombic b easy-axis, where the field-induced mode becomes Raman active (see 
discussion in the text). We observed a monotonic hardening of the field-induced mode with increasing applied field. 

Next we present the one-magnon Raman response in La2Cu04 for the (RR) polarization configuration. In this 
geometry, the electric-field of both the incident and scattered light is circularly polarized rotating anti-clockwise, 
Ef = Eg = (x Q + ix;,)/v / 2. According to the theory described in Ref. 0, the (RR) in-plane polarization configuration 
is the adequate Raman geometry in order to observe the field-induced mode for a magnetic field applied along the 
orthorhombic b easy- axis, because it probes directly the nonvanishing elements of the rotated xf Y matrix 



where 9 is the angle of rotation of the spin-quantization basis within the (be) planed The (RR) polarization configu- 
ration then probes the element i(F + H) sin 9 cos 9. The results are presented in Fig. Eland we observed a monotonic 
hardening of the gap with increasing magnetic field. Only the data points for higher fields are shown in Fig. 03 because 
the intensity of the peak rapidly drops for fields lower than 4 T. 

Finally we present the one-magnon Raman response in Sr2Cu02Cl2 also for the (RL) polarization configuration. 
The results are presented in Fig. 03 First, we observe that at zero applied field the Raman spectrum is continuously 
increasing up to the lowest accessible frequency of 2 cm -1 . As a consequence, we can only establish an upper bound 
of 2 cm -1 for the in-plane magnon gap. This very small value is consistent with the general expectation that the 
in-plane gap for Sr2Cu02Cl2 has a purely quantum origin. When the magnetic field is in the plane (main panel) 
one observes a hardening of the one-magnon peak, which allows us to clearly identify the magnon gap as the field 
increases, while no changes are observed for a magnetic field parallel to c (bottom inset). Observe also (top inset) that 
the signal corresponding to the two-magnon continuum (which starts at the edge of approximately twice the magnon 
gap) has much lower intensity with respect to the one-magnon peak. 



For a field applied parallel to the a orthorhombic axis, we obtain the conventional field dependence of the magnon 
gaps in a transverse fieldi the hardening of the mode in the field direction, while the second mode remains unchanged 



Here we indicate with m a ,m c the gaps of the magnon modes at zero field, cj o (0) = m a and uj c (0) — m c . Moreover, 
with respect to the preceding article?! we restored the physical units for the magnetic field, introducing the Bohr 
magneton /is = 0.4668 cm _1 T~ 1 and the gyromagnetic ratio for the field along a. We then obtain for the u> a mode 
the relation 




(15) 



V. MAGNETIC SPECTRUM IN La 2 Cu0 4 



A. H parallel to a 



u a (H) = y/ml + (gS/JLBH) 2 , 



LUc(H) = m c . 



(16) 



^1 = ml + la H 2 , la = ( 5 >b) 2 . 



(17) 



9 



T=10K (RL) polarization 




Raman shift (cm -1 ) 



FIG. 6: (Color online) One-magnon Raman response for SnCuOaCh at T — 10 K in the {RL) polarization configuration, where 
only the in-plane gap is Raman active. Main figure: hardening of the in-plane gap with increasing in-plane magnetic field. 
Insets: (top) zoomed immage of the main plot, showing the much weaker intensity of the two-magnons scattering, starting at 
approximately twice the magnon gap; (bottom) featureless response for perpendicular magnetic field (the out-of-plane (XY) 
mode is not Raman active in the backscattering geometry used here). 



As one can see in Fig. the data for w a (H || a) follow perfectly the previous equation, with a coefficient 7^ = 0.93 
(cm T) -2 estimated by fitting the data with Eq. ifTTIl . This value allows us to estimate the gyromagnetic ratio as: 

9s = yfyi/ftB = 2.0, (18) 

in good agreement with the estimate given usually in the literature. 
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FIG. 7: (Color online) Comparison between the experimental data and theoretical predictions for the field dependence of the 
DM gap for H || a. The closed circles are the experimental data. The line is the best fit obtained using Eq. flTfr. with -vf = 0.93 
(cm T)~ 2 . 
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B. H parallel to c 



When the field is applied along the c direction one observes a spin-flop of the ferromagnetic spin components along 
c, which are ordered antiferromagnetically in neighboring planes at low fiel d 1 ' 11 !?! 11 . This weak ferromagnetic (WF) 
transition has been indeed measured in Ref. (23 for the x — 0.01 doped compound, occuring at a temperature- 
dependent critical field of about 4 T at T = 50 K. Since the magnitude of the ferromagnetic spin components along 
c decreases with the temperature proportionally to the AF order parameter cr , also the critical field H c for the 
WF transition decreases with temperature. By properly taking into account the effect of quantum and thermal 
fluctuations on the order parameter one can determine the H C (T) curve reported in Fig. [Hl^ where we also sketch the 
spin configuration above and below the transition. However, as it has been discussed in Ref. Q], a first estimate for 
the critical field at low temperature and for the field dependence of the magnon gaps can be obtained by neglecting 
the renormalization of (Tq, Using the notation of Ref. 01 we will decompose the spin at site i in its staggered n 
and uniform L component, so that Si/S = e zC * Vi rii + aL. In the AF ordered state (n) = <jq, where in general cfq 
is corrected both by quantum and thermal corrections. At low temperatures, and neglecting quantum fluctuations 
((T = 1), one obtains that the critical field is 

where r] — 2JJ± is the energy scale associated to the interlayer coupling J±, and D + is the modulus of the DM vector, 
which controls also the H = value of the DM gap, D + = m^i. According to the analysis of Ref. Q, the magnon 




T(K) 



FIG. 8: (Color online) Phase diagram of the weak-ferromagnetic transition and spin configuration for H \\ c. The arrows with 
open tip represent the spins in each layer, the small ones on the right of the layers the uniform spin components. At low 
temperature and zero or low field the uniform components of the spins in each layers are ordered antiferromagnetically in the c 
direction. Above the critical field 11911 the uniform components order ferromagnetically in the c direction. As the temperature 
increases the uniform components decrease follwing the decrease of the AF spin component along 6, and as a consequence also 
the critical field for the WF transition decreases. The phase diagram has been calculated in Ref. Q (see Fig. 10) by using 
parameter values appropriate for La2CuC>4. 

gaps evolve, below H c , as 

w\ = m 2 a + 2g- sJ^ + {g c sl x B HD + f, 

"1 = m 2 c + (9>bH) 2 + 2 V - y/4 v * + (g°n B HD + ) 2 , H < H c , (20) 
while above the WF transition they are given by 



^ = ml + \g£nBHD + \, 

u? c = m 2 c + (g c sf i B H) 2 + \g c s fi B HD + \, H > H c . 



(21) 
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Using this set of equations we recognize that the parameter r\ can be determined from the jump of the gap at H c , 
since 



We can then estimate 



from which it follows also that 



jI(H-) =m 2 a - 2 V (V2 - 1), u 2 a (H+) = m 2 a + 2 V . 



V2 



ri = 2JJ± => j ± /j ~ 3.2 x kt 6 , 



(22) 



(23) 



where we used J = 130 meV. We can then recognize that the field evolution of the ui a gap above and below the WF 
transition is controlled by a single parameter 7 C 



uj 2 a = m 2 a + 2 V ~ VV + (7ci0 2 , H<H C 
= ml+lcH, H>H C , 

where 7 C is uniquely determined by 2r], H c and the gyromagnetic ratio of the c direction g c s 

2r) 



lc = 9 C sHbD+ 



(24) 



(25) 



As a consequence, we can extract g c s both by fitting the experimental data with Eqs. 112411 and obtaining a value 7^, 
and by using directly the last two equalities of Eq. In the former case we obtain 7^ = 21.83 cm~ 2 T~ 1 , which 

corresponds to g c s — 2.7. In the latter case instead we can use the value of 2ri obtained from the jump of the gap at 
H c , the value D + = m a = 17.48 cm -1 , and the measured value H c ~ 6.6 T, finding 



2 V 



2.65. 



Thus, both estimates give a value which is quite close to the one commonly quoted in the literature, g c s = 2.45. 



(26) 
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FIG. 9: (Color online) Comparison between the experimental data and the theoretical predictions for the field dependence 
of the DM gap for H || c. The closed circles are the experimental data. Observe that near the WF transition two values of 
the gap have been measured (see also panel (c) of Fig. 0J. Indeed, since the transition is first order in a finite sample the 
transition is observed as a coexistence of the two phases (above and below the transition) for a range of field values around the 
thermodynamic critical field. The lines are the fit using Eqs. {2U, with 7 f = 21.83 cm 1 



Observe that this value of g c s is only an upper bound because we did not include in the set of Eqs. il24"jl and in the 
definition of the critical field i|19H the quantum and thermal correction of the order parameter ctq, which reduce its 
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value with respect to cto = 1 used so far. As it has been explained in Ref. Q], when this effect is taken into account 
one must replace D + by D + /<r in Eqs. ll2fil) - lj2"T)l 



m„ 



2i] — \ Arf 



(g c sf i B H) 2 + 2r 1 -\U^ 



g c s H B HD^ 
0o 



H <H r 



(27) 



and 



g c s H B HD^ 



^ = m 2 c + {g c s fi B H) 2 + 



g c s fi B HD^ 



At the same time the critical field is a function of ctq according to 



He. 



Arj 



00 



g c s ^ B D+ 1 + a 2 



(28) 



(29) 



The order parameter oq is determined at each temperature by computing the effect of transverse spin- wave fluctuations, 
which depends on the magnon gap according to a general equation like 



o$(H,T) = l-I ± (u a (H),w c (H),T), 



(30) 



where the precise expression for I± is given in Ref. 01 Note that since the DM interaction introduces an explicit 
dependence of the magnon gaps on the order parameter, see Eqs. I|27ll - I|28|l . Eq. ll-'iOII is a self-consistency equation for 
(To. As far as the previous estimates of r\ and g c s are concerned, one can see that the solution of the full set of Eqs. 
ll27l - l|28l) f2H) and ij3f)l) can slightly modify the values previously obtained in two respects. First, the jump of the gap 
at the critical field will depend not only on r\ but also on (T - Second, the analogous of Eq. (l26l) will be 



4// 



00 



fi B H c D + 1 + erg 



(31) 



Observe that both these corrections will contribute to a reduction of g° with respect to the previous estimate. However, 
since the experimental determination of a in this sample is not available and the theoretical value depends on the 
approximations used to derive Eq. 1(3(1)1 , we retain here the estimates given above using cto = 1 and we refer the reader 
to Ref. Q] for a more detailed discussion of this issue. 

We should point out that a somewhat larger estimate for the critical field H c has been extracted recently from the 
neutron scattering measurements of Ref. [13, in a La2Cu04 sample with almost the same Neel temperature as the 



one considered here. According to the previous Eq. ll29f) . several factors can affect the critical field. Thus, the larger 
critical field measured in the sample of Ref. can be explained with a larger value of the staggered order parameter 
and of the interlayer coupling rj, or also with a smaller value of the DM vector D + . However, one should remember 
that the weak-ferromagnetic transition for H \\ c is a first-order transition accompanied by a hysteresis that can be 
large at low temperature (see for example Ref. 113), and hence the experimentally determined critical field could be 
affected by the hysteresis. As a consequence, further investigation of the parameter values for this sample can shed 
more light on this apparent discrepancy between the Raman and the neutron scattering measurements. 



C. H parallel to b 

Finally, let us discuss the case of H || b, where also the data for the field-induced mode u c are available for H > 4 
T. Here again the field dependence of the magnon gaps follows a different behavior for small and large field. As it 
has been explained in Ref. Q, due to the DM interaction a field along b generates an effective staggered field in the c 
direction, giving rise to a rotation of the staggered order parameter in the be planed As a consequence, at low field 
the classical configuration of the staggered AF order parameter n m in the m-th plane is given by 



n m = (0,(T b ,(-l) m (7 c ). 



(32) 
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where <j\, and a c indicate, respectively, the components of the order parameter along b and c direction. Observe that an 
oscillating staggered <r c component implies that in the spin decomposition Si/S = e lCJ ' ri n i + Lj given above actually 
the S c components coming from n c are ordered ferromagnetically in neighboring planes, to allow for the (average) 
uniform spin components L induced by the DM vector to align along the applied fieldi 

As the field strenght increases, the spins perform a two-step spin-flop transitionti^ii at H — the in-plane AF 
components allign along the a direction, while at H = > He the a component vanishes and only the staggered 
spin components in the c direction are left. In Fig.E3 we shown the (H — T) phase diagram which has been evaluated 
m Ref. |lj by means of a saddle-point approximation for the transverse spin fluctuations. The temperature Tn is the 
one where the in-plane spin component (along b or a) vanishes, leaving a non-zero <r c component. Observe that the 

jump of T N {H) at HP is an artifact of the approximation: indeed, the critical field for this transition is controlled 
by the energy scale where the u> a gap vanishes (see Eq. 13411 below), i.e. 

g b ji B HP =m a = D +1 (33) 

which corresponds, in our case, to — 18 T. Since within the saddle-point approximation m a is temperature 

independent, at this level one cannot recover the temperature variation of He which would eliminate the anomalous 
jump of Tjv(-ff) reported in Fig.EB Since the maximum field used in the present experiments is H = 11 T, this spin 
flop transition is out of the range accessible in our measurements. It is perhaps worth pointing out that in a recent 
work the evolution of the magnetic signal for H || b has been measured by means of neutron scattering up to fields 
of 14 Tiii By assuming a power-law decay of the intensity of the (100) neutron peak - corresponding to magnetic 
order along the b direction - the authors could estimate for the above critical field a slightly higher value of 22 T. 
Considering that a theoretical prediction of the neutron-peak intensity as a function of magnetic field is still lacking, 
this estimate seems in good agreement with the value ll33|l given above. 
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FIG. 10: (Color online) Phase diagram and evolution of the spin configuration for H \\ b. For the spin components we use the 
same notation as in Fig. Elabove. At H < He both the spins and the uniform components rotate in the be plane, so that the 
net uniform magnetization in along the b axis. In this region Tm(H) is the temperature at which the in-plane AF component 
along b vanishes. At H^ < H < H^ the AF spin components allign along a which is perpendicular to the plane of the figure, 
so that these spin components, which vanish at the temperature Tn(H), are represented by the small crosses. At the same 

(2) 

time, the uniform components align ferromagnetically in neighboring planes. At H > He only the c components survive. 

As far as the magnon-gap evolution with magnetic field is concerned, at field below it is the usual one for 
longitudinal fields^ 

«l - + {9 b s , B Hf + ^(^ y ^) 2 + 4 (5 ^) 2 (5±^). (34) 
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At small field these expressions can be approximated as 

ul = ml- la H\ 7a = (g b sf i B ) 2 Am a , Am a 

uj 2 c = m 2 c + lc H 2 , 7c = (g b s n B ) 2 Am c , Am c 

Using Eq. l(34T) and the zero-field value of the a mode we obtain again an excellent agreement with the experimental 
points, as one can see in Fig. El For the sake of completeness we also report in Fig. El the approximate expression 
(13511 . which is indeed valid until the field Hi ~ 4 T. From the fit of the uj a data using Eqs. (13411 we obtain both 
m c = uj c (H = 0) and the gyromagnetic ratio. The results are 

m c = 36 cm -1 , g b = 2.08, (36) 

which are again in excellent agreement with the values reported in the literature. 

In the inset of Fig. El we show also the field dependence of the field-induced uj c (XY) mode according to Eq. Q34H , 
using the values l.Uill extracted from the fitting of the uj a mode. As one can see, only the data point at 9 T is out of 
the range of the thoeretical curve. The appearence of the field-induced ui c mode for the in-plane Raman scattering 
in the (RR) polarization configuration is now well understood as being a consequence of a continuous rotation of the 
magnetization axis when H || 5il& 
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FIG. 11: (Color online) Comparison between the experimental data and the theoretical predictions for the field dependence of 
the a mode for H || b. The blue solid line is the fit done using the first of Eqs. 13411 with g s and m c as fitting parameters, giving 
g b a — 2.08 and m c = 36 cm -1 . The green dashed line is the approximated expression 13"it evaluated with the same parameters 
values, and it is only valid at low field. Inset: field dependence of the c mode observed at fields larger than 4 T. The solid line 
is the curve corresponding to the second of Eqs. (JSU- 



D. Peak intensity as a function of magnetic field 

A complete microscopic theory of the Raman scattering which allows one to compute the exact shape and the 
absolute intensity of the Raman peaks is out of the scope of the present paper. Nonetheless, when data taken at 
different magnetic fields are compared, one could expect that the main dependence of the relative peak intensity 
measured by Raman can be at least qualitatively described by the theory developed in Ref. 01 Indeed, the spectral 
function A a , c contains itself an intrinsic dependence of the peak intensity on the magnetic field, which enters in the 
Raman response trough the relation 11311 . Thus, we can evaluate Eq. 113t using the theoretical prediction for A a , c and 
compare it with the experimental data. As we shall see, even though our analysis does not include the spin damping 
processes, the overall agreement between the theoretical predictions and the experimental data is fairly good. 



ml + mi 



(35) 
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As it has been discussed in Ref. Q], the spectral function of each mode is defined from the Green's function for the 
corresponding fluctuations. In the absence of magnetic field the Green's function for the transverse mode is a diagonal 
matrix 



G 1 - ~ ( n aK ' ,,,2x^2^ ), (37) 



i (rt + 4(k) o 



gc 



where the magnon gap is by definition iv a _ c — e a ,c(k = 0). In this case the Green's function matrix is also diagonal 
and we simply obtain 

A a (u > 0) = -~ImG a {iu) n -> ui + i0+) = ^—S(uj - u> a ), (38) 

TT ZLO a 

and analogously for A c . Thus, from Eq. I|38ll we can deduce that the peak intensity I a evolves as l/u> a , i.e. it is larger 
for smaller gap values. 

When a magnetic field is applied, two different cases must be considered: (i) if the matrix G is still diagonal the 
structure (138)) of the spectral function is preserved, and both the peak position and its intensity evolve as uj a , c (H), 
and I a = l/u> a ,c(H), respectively; (ii) if off-diagonal terms proportional to the magnetic field appear in Eq. i|37H the 
Green's function of the transverse fluctuations has a non-diagonal structure which leads to the appearance of several 
magnon gaps in the response of each single mode. For example, in the case of longitudinal field one has for the spectral 
function of the a mode a structure likei 

A a {cj > 0) 



ZUJ n ZLU r 



(39) 



However, one finds in general that Z a /uj a ^> Z c /oj c , so that essentially a single peak at the uj a energy is observed in 
the measurements, but the spectral weight of this peak is I a = Z a /uj a and not just l/u> a as in Eq. I|38|) . Indeed, the 
factor Z a in Eq. 13911 leads in general to an additional field dependence of the intensity on the magnetic field. 

For an ordinary easy-axis AF Eq. l(38jl is valid when a transverse field is applied, so that an hardening of the gap in 
the field direction should be also accompanied by a softening of the peak intensity. For a longitudinal field one finds 
instead the spectral function H39J1 given above. According to the calculations of Ref. Q], Z a is given by 

-u? n + ml - H 2 , s 

Z a = 2 2 , 40 

and it is a decreasing function of the magnetic field. As a consequence, the overall factor Z a /uj a in Eq. (139)) is 
increasing much more slowly than the ~ l/oj a behavior that one could expect for a transverse gap, as it is shown in 
Fig- El Here, according to Eq. 11311 . we included also the bose factor [ns(w) + 1], whose contribution to the overall 
dependence of the peak intensity is however very small. 

In the case of H || c in an ordinary easy-axis AF one would expect the peak for the a to be unchanged, as indeed 
observed in S^CuCbC^. However, for La2CuC>4 the presence of the DM interaction does modify both the peak 
position and its intensity as a function of magnetic field. Below the critical field for the spin-flop transition the 
spectral function of the a mode ial 

A.(«>0)w— 5{u-w a ), (41) 



where 



Once again, while the l/w a (H) increases with magnetic field the Z a factor decreases, giving rise to an almost constant 
spectral weight below the spin-flop transition, see Fig.^3 Here we used the r\ value extracted from the fit with <jq = 1. 
Above the spin-flop transition the matrix of the transverse fluctuations is again diagonal and the standard I a = 1/ ui a 
spectral weight is expected, with the uj a (H) dependence given by Eq. l|28jl . 



VI. MAGNETIC SPECTRUM IN Sr 2 Cu0 2 Cl 2 



The case of Sr 2 Cu02Cl2 is rather simple. This system is very well described by the Hamiltonian Q with a large 
XY anisotropy a c ~ 10~ 4 J. Since the system is tetragonal, one does not have in principle any intrinsic in-plane 
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FIG. 12: (Color online) Field dependence of the normalized intensity I a (H)/I a of the a peak for a longitudinal field. We 
estimated the relative error on the measured peak intensity around 20 %. Solid line: field dependence of the peak intensity 
according to Eqs. 1.3t)ll and 14011 . Observe that by neglecting the contribution of the Z a factor from Eq. 14011 one would obtain 
I a = 1/i^a, corresponding to the dashed line, which shows an increasing of the relative intensity not observed in the experiments. 



tr 0.6 




FIG. 13: (Color online) Field dependence of the intensity of the a peak for a transverse field in the c direction. Observe that 
the spin-flop WF transition is first order, so that by increasing the magnetic field one first access a crossover regime where the 
Raman intensity is split between the two peaks corresponding to the equilibrium gap values at H < H c and H > H c . Thus, 
only the data at H > 7 T should be compared with the dashed curve, corresponding to the state above the spin flop. 



anisotropy, i.e. one would expect a a = in this case. Nonetheless, it has been suggested in Ref. @ that a very 
small in-plane gap can still exist as due to purely quantum effects. Recent ESR measurements seem to confirm this 
prediction and give an estimate of the in-plane gap as m a = 0.048 <~ 0.05 meV.- As we can see from Fig. there 
might be a peak at low energies (below 2 cm" 1 for H = T), which can well correspond to the same in-plane gap as 
observed by ESR. As far as the field dependence of the gaps is concerned, the results obtained in the Sec. III-IV of 
the preceding article for a conventional easy-axis AF still apply,- Moreover, since no orthormonbic distortion exists in 
Sr 2 Cu02Cl2, no DM interaction is present, and no effects due to staggered fields will show up in the field dependence 
of the magnon gaps. 

In Fig. El we compare the extracted dispersion of the in-plane gap in Sr2Cu02Cl2 for an applied magnetic field 
parallel to the CuC>2 layers, with the predictions for a transverse magnetic field from the theory presented in the 
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preceding article^ 

u a (H) = y/ml + (g st i B H)\ (43) 

As we can see, by using the value of the gap given by the ESR measurements, m a — 0.048 ~ 0.05 meV, 6 and g s = 2.05 
the agreement is already quite good. On the other hand, we can use the Raman data alone to give an independent 
estimate of both m a and g s , in the same spirit of the previous Sections. This is shown in Fig. EI and the best fit is for 
g s = 1.98 and m a — 1.96 cm" 1 (m a — 0.24 meV). Observe that this estimate of m a is larger than the one measured 
by ESR, but it is still consistent with the fact that at H = no gap was observed in the Raman spectra at frequency 
larger than 2 cm" 1 , which is the lower bound of the accessible frequency range as shown in Fig. El 
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FIG. 14: (Color online) Field dependence of the in-plane spin-wave gap measured by Raman in Sr2Cu02Cl2. Solid line: field 
dependence of the a gap obtained using the value m a = 0.4 cm - and g B = 2.05 measured by ESR and the theoretical prediction 
for a transverse field. Dashed line: best fit on the data point using Eq. 14311 with m a and g s as fitting parameters. 

Finally, we should emphasize that no field-induced modes were observed in Sr 2 Cu02Cl2 neither for perpendicular 
nor in-plane magnetic fields. This results from the absence of the Dzyaloshinskii-Moriya interaction in Sr2Cu02Cl2- 
In fact, it is exactly the DM interaction in La2Cu04 that causes the rotation of the spin-quantization basis and 
modifies the Raman selection rules allowing for the appearance of a field-induced mode when the field is along the 
easy axisiA& Moreover, the DM interaction is responsible for the in-plane gap evolution for a field along c, which is 
instead unexpected in a conventional easy-axis antiferromagnet. 

VII. CONCLUSIONS 

We have measured, discussed, and compared the one-magnon Raman spectrum of Sr2Cu02Cl2 and La2Cu04. We 
have seen that, for the case of Sr 2 Cu02Cl2, which is a conventional easy-axis antiferromagnet, there is an in-plane 
magnon mode at very low energies in the DA g channel that is accessible in the (RL) polarization configuration. No 
out-of-plane nor field-induced modes were observed in the geometries used, in agreement with the general expectation 
for a conventional easy-axis antiferromagnet. For the case of La 2 Cu04, the magnetic field evolution of the in-plane gap 
measured in the (RL) polarization configuration is made rather nontrivial due to the presence of the DM interaction. 
When the field is along c the DM gap first softens, jumps discontinuously at the critical field for the WF transition and 
then hardens. At the same time, when a longitudinal field is applied the a mode softens and the (field-induced) c mode 
appears in the (RR) polarization, as a consequence of a rotation of the spin quantization basisA The magnetic-field 
dependence of the one-magnon Raman energies were found to agree remarkably well with the theoretical predictions 
of part I of this workji allowing us to extract from the Raman spectra the values of the various components of the 
gyromagnetic tensor and of the interlayer coupling. Moreover the analysis of the field-evolution of the Raman-peak 
intensity, which also shows a general good agreement with the experiments, demonstrated that the long- wavelength 
analysis of Ref. 1 contains new useful informations with respect to the standard spin- wave calculations of the magnon 
gaps existing in the literature. 
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